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Abstract 

The electric polarizability and logarithmic mean-excitation energy are calcu- 
^ . lated for the deuteron using techniques introduced in atomic physics. These 

rH I results are then used to improve limits on the atomic-deuterium frequency shift 

due to nuclear polarization in the unretarded dipole limit, as well as conhrming 

previous results. 



The remarkable experiments [|I], 0] currently being performed on the isotope shift in 
atomic hydrogen (^H vs. ^H) are primarily determined by differences in the masses of 
the isotopes, but are significantly sensitive to nuclear structure. These measurements 
provide the most precise determination of the difference in sizes of these isotopes. 
The most recent result for the d-p isotope shift in the 1S-2S level splitting is 

Az/D_H = 670994334(2) kHz , (1) 

of which roughly 5000 kHz is attributable to the finite-size differences of the nuclei, 
while roughly 20 kHz is due to the electric polarizability of the deuteron. In other 
words, in addition to a weaker Coulomb potential arising from the nuclear charge 
distribution seen by the electron at very short (on the atomic scale) distances, that 
electron also "distorts" or polarizes the nucleus, which enhances the binding. Four 
numerical calculations of the effect of nuclear polarization on the isotope shift have 
been performed recently p, |, |^, |^, although the relevant leading-order analytic 
results for the nth S-state have long been known [0, H: 

AEpoi = -5me a |0„(O)p «s (^ + log (^) ) ' (2) 

where a is the fine-structure constant, me is the electron mass, |0„(O)p = fi^a^/irn^ 
is the square of the wave function of the electron at the origin, oe is the deuteron 
electric polarizability, /i is the e — d reduced mass, and we work in natural units 
{h = c = 1). Even though uncertainties in the polarization calculations are currently 
smaller than the error quoted in Eq. (1), planned improvements^ in that accuracy 
warrant a strong effort to reduce the theoretical uncertainty to a minimum. 



The electric polarizability of a nucleus (or atom) is defined by |T(]] 



".-fi:™^. (3) 

where Eq is the energy of the ground-state |0), E^- is the energy of the Nth excited 
state, and D is the electric-dipole operator, which effects the transition between those 



states. The definition (3) can be rearranged into the form of a sum rule [p!Ot ^ 
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where a!^'^{uj) is the cross section for photoabsorption of unretarded-dipole (long- 
wavelength) photons by the nucleus. Concomitantly, the logarithmic mean-excitation 
energy in Eq. (2), E, is defined by 

2a ^ \{N\D\^ j^^ ^^^^ _ ^^^^^j ^ ^^ log(E/m,) , (5) 

and clearly corresponds to placing a factor of log(c(j) in the integrand in Eq. (4). 
The august o"_2 sum rule and its (less well-known [O]) logarithmic relative aig have 



been used to evaluate aE and E by explicitly constructing {N\D\Q) (or equivalently, 
CT^'^{uj)) and performing the integral numerically. Results for ue for many "realistic" 
potential models are known, although the two most recent calculations @, H did not 
have any models in common. 

In this work we will: (1) calculate ue for a set of models that subsumes most of 
those of Refs. [Q and P and includes several more; (2) calculate log [E) for these 
models; (3) use novel (for nuclear, but not atomic, applications) numerical techniques 
for calculating both aE and log (-E); (4) critically discuss the potential models and 
attempt to assign a subjective but credible uncertainty to the results. In this way 
we will confirm the previous results, while shrinking the uncertainty associated with 
them. Our numerical techniques were first applied to atomic problems, but now find 
a comfortable home in nuclear physics. 

The technique we use for calculating ue was first used by Podolsky [0 to treat 
dispersion in hydrogen atoms. Definition (3) is fully equivalent to 

as = 2a(0|D,|A^,), (6) 

where 

(if-Eo)|A^,) = DM (7) 

is solved subject to finite boundary conditions. Note that D does not connect the 
ground state (the only bound state) of the deuteron to itself. Resolution of Eq. (7) 
into partial waves, incorporation of the nuclear (including tensor) force, and other 
minor (although tedious and important) details are contained in Ref. ||I0|, together 
with many analytic results for simple potentials. Our calculation will employ the 
usual nonrelativistic (classical) dipole operator. 

The resulting procedure is only slightly more complex than solving for the deuteron 
ground state, and it is very stable. We have calculated aE for 14 different "realistic" 
nucleon-nucleon (NN) potential models. Such models must contain OPEP (the one- 
pion-exchange potential), which dominates the binding of light nuclei, and they must 
fit the NN data reasonably well. All of the models used in Refs. 0] and are in 
this category, although the quality of the fits (of various potential parameters) to 
the data differs rather dramatically from case to case. Most of those models could 
be characterized as "first-generation" models. Recently, the Nijmegen group and 
their collaborators |T^ , p^ , |T6[| have constructed "second-generation" models, which 



provide good- to very-good-quality fits to all NN data, even approaching x^ per degree 
of freedom ~ 1. Such fits are sufficiently good that they can be regarded as alternative 
phase-shift analyses. This does not necessarily imply that the underlying physics has 
a corresponding accuracy, since several of these models are purely phenomenological, 
except for the all-important and dominant OPEP that incorporates different pion 
masses in different (isospin) states. 

We determine the logarithmic mean-excitation energy E (or logarithmic sum rule) 
using a trick developed for calculating various logarithmic mean-excitation energies 
in atoms |]T^, one of which is the Bethe logarithm. If we add a parameter X = C, ■ f 
to {H — Eq) in Eq. (7), where ^ is dimensionless and / has the dimensions of energy. 



we can then define (and easily calculate) 

(^E[i) = — 2^ 7T— ^ ^, (8) 

where aE^O) is the usual result. The integral of ctEiO fro^i to A (very large 
compared to E/f) generates 

f\cMO oc - E '^^'^'°^'' log [(En - Eo)/Af] , (9) 

which gives the desired logarithm. A similar integration gives 

r|„,«)oci:^^iog[(^.-i.o)A/i. (10) 

For numerical purposes, we split the integral in Eq. (10) into /^ +/i°°, and the di- 
mensionful scale parameter / determines where the split occurs in energy units. Re- 
arranging slightly and changing variables to 1/^ in the second integral, we achieve 
our final result 

as(0) log {2E/m,) = /' ^[aeiO ' «£;(0) + «e(1/0] - «i?(0) log(m,/2/) . (11) 
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The integrand is finite everywhere. Choosing / ~ 3 — 5 I-EqI rnakes the integral 
converge to 5 significant figures with only a few (~ 6) Gauss quadrature points, and 
all results for E are independent of / if the integrals are performed with sufficient 
accuracy. Podolsky's method [l^, |T^ makes the calculation of q;_b(0 as easy as that 
of Q!e(0). The method is very stable. 

Table 1 presents our results for ue, \og{2E/me), and Ai^poi separated into first- 
generation |]18|, |19|, ^, |21|, ^, ^ |2^ (listed in order of appearance in Table 1) and 



second-generation [14, 15, |T6| (potential) categories. Note that there is much more 
spread in the first-generation results, reflecting indifferent fits to the NN data. The 
spread in the second-generation results can be summarized by 

z/poi = 19.26(6) kHz, (12) 

and 

aE = 0.6328(17) fm^ . (13) 

As noted below in Ref. 0], this is not a numerically complete result for the sum of 
all polarizability corrections, since it incorporates only unretarded dipole approxima- 
tion. Higher multipoles, retarded dipole contributions, seagulls, etc., have not been 
included here, and may decrease this result by up to 1 kHz [3. 

All of the appropriate results are quite close to those previously calculated [^ 0, 
with the electric polarizabilities differing at most by 2 in the last quoted significant fig- 
ure in those references. Such small differences could be attributed to slightly different 



Table 1: Deuteron electric polarizabilities, oe, in units of fin^, logarithmic mean- 
excitation-energy ratios, log(2£'/me), and deuteron 1S-2S polarization-energy shifts, 
z/poi, in kHz. 



Potential Model 


aE(fm^) 


\og{2E/m,) 


z/poi(kHz) 


Second-Generation Potentials 


Reid Soft Core (93) 


0.6345 


2.9616 


19.31 


Argonne Vis 


0.6343 


2.9625 


19.31 


Nijmegen (loc-rel) 


0.6334 


2.9618 


19.28 


Nijmegen (loc-nr) 


0.6327 


2.9624 


19.26 


Nijmegen (nl-rel) 


0.6328 


2.9619 


19.26 


Nijmegen (nl-nr) 


0.6319 


2.9625 


19.24 


Nijmegen (fuU-rel) 


0.6311 


2.9615 


19.21 


First-Gene 


ration Potentials 




Reid Soft Core (68) 


0.6237 


2.9638 


18.99 


Bonn (CS) 


0.6336 


2.9630 


19.29 


Paris 


0.6352 


2.9627 


19.34 


de Tourreil-Rouben-Sprung 


0.6376 


2.9623 


19.41 


Argonne V14 


0.6419 


2.9624 


19.54 


Nijmegen (78) 


0.6472 


2.9612 


19.70 


Super Soft Core (C) 


0.6497 


2.9617 


19.77 



versions of the potentials (new potentials are often a matter of "work in progress'' 



Note also that the pairs of new Nijmegen |15] local and nonlocal potentials (labeled 
"loc" and "n/" in Table 1) have versions with relativistic ("rel") and nonrelativistic 

2M 



.2 
■rer 



("nr") kinematics (corresponding to identical deuteron energies of IJM"^ — k^ 

2 

or — ^, respectively). The slightly smaller value of Krei in the (excellent) zero-range 
approximation [^ |T^ accounts for those differences in the values of aE, although this 
makes relatively little difference in i/poi. The "full" Nijmegen potential |^ has the 
same form in all partial-waves and fits the NN data less well than the others. 

The result (13) agrees very well with a prediction [|1^] of aE = 0.632(3) fm^ made 
many years ago, and this warrants further comment. One can perform perturbation 
theory about the "zero-range" limit by turning off the forces in p-waves, dropping 
the deuteron d-state, and replacing the (reduced) deuteron s-state wave function 
by its asymptotic form: u{r) = Asexp{—Kr), where As is the s-wave asymptotic 
normalization constant. With this ansatz we obtain [^ |lOl 



aE 



~ ^0 

= Or 



afiA's 
32fi;5 



(14) 



where fi is here the n-p reduced mass, and log{2E/'mf.) = 2.9671. This remarkably 
simple formula overestimates the complete result by approximately 1% . There is little 
uncertainty in any of the quantities except for As, which was recently determined [p5 



to be As = 0.8845(8) fm"^'^ in agreement with the value used in Ref. |Ty], and which 
leads to a% = 0.6378(12) fm^. Moreover, the corrections, Aue, to a% defined by 
«£ = «£; + Aa^; can be determined from the potential models (see Refs. p and |1T0(| ) 
to be AaE = -0.0044(2) fm^, which leads directly to ue = 0.6334(14) fm^, which is 
consistent with Eq. (13). Note that no relativistic corrections have been incorporated 
and they are not likely to be negligible on the scale of the uncertainty in Eq. (13). 

Why do the "second-generation" potentials agree so well with the perturbation 
theory estimates? The answer is that As is determined by analyzing NN scattering, 
and we stated earlier that the new potentials could be viewed as alternative phase- 
shift analyses. That is, they fit the NN data quite well, and associated properties 
(such as As) should agree with other experimental determinations. Thus, aE is very 
well determined. 

We summarize by noting that oe and log(-E') have been calculated for the deuteron 
by novel methods. These calculations confirm previous results and add additional 
ones. We strongly recommend that only second-generation potential results be used 
when assessing the reliability of oe calculations. Equation (12) gives our best esti- 
mate for the leading-order (unretarded-dipole or long-wavelength) approximation to 
the nuclear-polarizability correction given by Eq. (2). 
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